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In this paper we study for a given azimuthal quantum number k the eigenvalues of the 
Chandrasekhar-Page angular equation with respect to the parameters /i := am and v := au, where 
a is the angular momentum per unit mass of a black hole, m is the rest mass of the Dirac particle and 
to is the energy of the particle (as measured at infinity). For this purpose, a self-adjoint holomorphic 
operator family A(n; fi, v) associated to this eigenvalue problem is considered. At first we prove 
that for fixed « 6 M \ (— J, |) the spectrum of A(k\ fj,, v) is discrete and that its eigenvalues depend 
analytically on ([i,v) £ C 2 . Moreover, it will be shown that the eigenvalues satisfy a first order 
partial differential equation with respect to /i and v, whose characteristic equations can be reduced 
to a Painleve III equation. In addition, we derive a power series expansion for the eigenvalues in 
terms of v — fi and v + fj,, and we give a recurrence relation for their coefficients. Further, it will 
be proved that for fixed (fi, u) £ C 2 the eigenvalues of A(k; h, v) are the zeros of a holomorphic 
function O which is defined by a relatively simple limit formula. Finally, we discuss the problem if 
there exists a closed expression for the eigenvalues of the Chandrasekhar-Page angular equation. 



I. INTRODUCTION 



The angular eigenvalue problem of a spin-i particle in the Kerr-Newman geometry is given by the Chandrasekhar- 
Page angular equation (see 0, Chap. 10, Sec. 104]) 

CtS + i=(amcos9-X)S_L, (1) 

2 2 2 

CTS_i = (am cos + A) (2) 

where the Kerr parameter a is the angular momentum per unit mass of a black hole and m is the rest mass of the 
Dirac particle. Moreover, the differential operators £f are defined by 

2 

2 2 sin 

where u> is the energy of the particle (as measured at infinity) and k is a half- integer, i.e., n = k— i with some k € Z. 
A parameter A € R is called an eigenvalue of this spectral problem if the system given by - © na s a nontrivial 
solution which is square-integrable on (0, it) with respect to the weight function sin (9. In this paper we study for fixed 
k the eigenvalues of the Chandrasekhar-Page angular equation as a function of the parameters fi := am and v := au. 
As a main result, we will prove that the eigenvalues satisfy a first order quasi- linear partial differential equation, and 
we will derive a power series expansion for the eigenvalues in terms of v — fi and v + fx. 

For this purpose it is necessary to consider the system Q - J2J in a more general context where k is real, \k\ > 5, 
and /j, v are complex numbers. At first we rewrite this system for fixed k£I\ 5) as an eigenvalue problem for 
some self-adjoint holomorphic operator family A — A(k; \i, v) depending on the parameters v) € C 2 . In the special 
case where (p, v) G M 2 the differential operator A(k; /1, v) is self-adjoint and has purely discrete spectrum. In Section 
II we prove that for a given k the eigenvalues Xj (k; fj,, v) of A are holomorphic functions in v), and we derive some 
basic estimates for them. Furthermore, we transform the system Q - (J5J) to a matrix differential equation 



y'{x) 



-B Q + — !— B 1 + C 

x x — 1 



!/(*) (3) 
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on the interval (0, 1) with coefficient matrices 

V o i + i) -S-j/ 

which can be extended to the complex domain C \ {0, 1}. In this way we obtain a further characterisation of the 
eigenvalues of A and some useful estimates for the corresponding eigenfunctions. Applying analytic perturbation 
theory, we show in Section III that the eigenvalues Xj(n; fi, v) satisfy the partial differential equation 

<9A <9A 

{u-2vX) — + {v-2u\) — +2kll + 2uv = Q. (4) 

Of! OV 

In particular, this result can be used to obtain a recurrence relation for the coefficients c m ^ n of a power series expansion 

oo 

Xj(n;fi, v)= ^2 c m.n{. v - y) m { v + A*)™- 

m,n— 

In Section IV we solve the PDE J2J by the method of characteristics. First, we derive an explicit formula for the 
eigenvalues in the case \fi\ — Moreover, in the regions where ^ \v\ we reduce the characteristic equations of 
Q to a Painleve III equation 

vv' + t v v" - t(v') 2 - 2 k (v 2 ± l) v - t (V 4 - l) = 

with parameters a — ±/3 = 2k and 7 = — 5 = 1 according to the notation in 0] and 0|. As this differential equation 
is in general not solvable in terms of elementary functions, we cannot expect a closed expression for the eigenvalues of 
the Chandrasekhar-Page angular equation for all (fj,, v) el 2 . However, if k is a half-integer, i.e., n = k — ^ with some 
positive integer fc, then a±(3 — 2{2k— 1), and there are integrals of polynomial type for the third Painleve equation in 
this special case, cf. 0]. Hence, if n — ±i, ±|, . . ., there exist algebraic solutions of the partial differential equation 
(@J, and the question arises if these explicit solutions are in fact eigenvalues of the Chandrasekhar-Page angular 
equation. It turns out that there is another type of "special values" associated to the operator A, called monodromy 
eigenvalues, which belong to the algebraic solutions of the PDE For a half- integer k, the monodromy eigenvalues 
are introduced in Section V by requiring that the system @ has a fundamental matrix of the form 

[x(l - x)}~%~i H(x) 

with an entire matrix function H : C — > M2(C). This property turns out to be equivalent to the existence of special 
solutions of the form 

[x{l-x)]~%-ip ± {x)e ±2tx 1 



where : C — ► C 2 are polynomials and t = ±^/v 2 — p? . For comparison purposes, an eigenvalue of A can be 
characterised by the property that J3J possesses a nontrivial solution of the form 

- x)\% + ^r){x) 

with some entire vector function r\ : C — ► C 2 . We prove that the monodromy eigenvalues are zeros of a polynomial 
with degree 2k — 1 whose coefficients are polynomials in /1 and v. Moreover, it can be shown that monodromy 
eigenvalues and "classical" eigenvalues are distinct at least in a neighbourhood of (fi, v) = (0,0). Nevertheless, they 
are both characterised by the fact that certain monodromy data of the system are preserved for all parameters 
(/i, v). In fact, A is a monodromy eigenvalue of A if and only if the monodromy matrices of 10 at the regular-singular 
points and 1 are diagonal, whereas A is a classical eigenvalue of A if and only if a certain non-diagonal entry of 
the connection matrix for the fundamental matrices at and 1 vanishes. Hence, for the Chandrasekhar-Page angular 
equation the monodromy as well as the classical eigenvalue problem is closely related to the isomonodromy problem 
for the differential equation Q- Monodromy preserving deformations for such a system were studied by Jimbo, Miwa 
& Ueno in however, only for the case that the eigenvalues of Bo and B\ do not differ by an integer, i.e., n + | ^ Z. 
In Section VI we consider the isomonodromy problem for JJJl in the case that k is a half-integer. As a consequence, we 
show that the monodromy eigenvalues of A satisfy the partial differential equation and we obtain an alternative 
derivation of (0J for the classical eigenvalues of A. Unlike the proof in Section III, which relies on the particular 
structure of the Chandrasekhar-Page angular equation, the method presented in Section V is more general and based 
on finding suitable deformation equations for parameter-dependent differential equations. Thus, we expect that this 
technique is applicable to other eigenvalue problems as well. 
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II. A SELF-ADJOINT HOLOMORPHIC OPERATOR FAMILY ASSOCIATED TO THE 
CHANDRASEKHAR-PAGE ANGULAR EQUATION 



By introducing the notations 



/i := am, v := aw, :~ V sin 

the Chandrasekhar-Page angular equation JQ) - (J2J) takes the form 
W*H«) := I , „ \S'(9)+( 




sin 9 



^COStf 

v sin 



-b — v sin f 

jJL COS 



S(0) = \S(6), 6>e(0,7r), 



(5) 



with fixed k € K \ (—5,3) and parameters (/x, z^) € C 2 . We can associate the so called minimal operator Aq to the 
formal differential expression 21, which acts in the Hubert space Ti := Jzf 2 ((0,7r),C 2 ) of square integrable vector 
functions with respect to the scalar product 



(Si, S2) 



S 2 {9)*S 1 (9)d9, S u S 2 eH. 



(6) 



The operator A given by V(A ) = Cg° ((0,tt),C 2 ) and A S := 215 for S £ T>(A Q ) is densely defined and closable. 
For I I > i and (/i, v) £ ]R 2 the formal differential operator in (JSJ) is in the limit point case at and 7r, hence A 
is even essentially self-adjoint. In the following we denote the closure of Aq by A = A(k,;^jl,v). According to 
Theorem 5.8] the domain of A(k; 0, 0) is given by 

V(A) = {S £ U : 5 is absolutely continuous and A(k; 0, 0)5 £ H} . 

Since A(k; /1, v) = A(k; 0, 0) + T(p, v) with the bounded multiplication operator 



T([i,v) 



-/icost' -i/sinf 
-vsm.9 )icos9 



its domain of definition T>(A) is independent of (fj,, v) £ C 2 (see jQ, Chap. IV, §1, Theorem 1.1]). Moreover, if 
(ji, v) el 2 , then T(jj,, v) is a symmetric perturbation of A(n; 0, 0), and 0, Chap. V, § 4, Theorem 4.10] yields that 
A(k;^l,v) is self-adjoint. Thus, according to Kato's classification [lj, Chap. VII, §3], A(k;^l,v) forms a self-adjoint 
holomorphic operator family of type (A) in the variables (/i, v) £ C 2 . Further, the spectrum of A(k;0,0) is discrete 
and consists of simple eigenvalues given by 



A J («;0 ) 0)= S gn(i) |k| 



Ul 



j G z \ {0} 



(7) 



(for the details we refer to Appendix A). This means, in particular, that A(k;0, 0) has compact resolvent, and from 
(I4L Chap. V, § 2, Theorem 2.4] it follows that A(k; \i, u) has compact resolvent for all (/x, v) £ C 2 . As a consequence, 
the spectrum of A(k; /z, v), (/j,, v) £ C 2 , is discrete, and since A(k; /j,, v) is in the limit point case at 9 — and 9 = tt, 
it consists of simple eigenvalues for (/1, v) £ R 2 . Now, [l4l Chap. V, §3, Theorem 3.9] implies that the eigenvalues 
Xj = Aj (k; /i, j/), j £ Z\{0}, of A(k; /j,, v) are simple and depend holomorphically on (/i, v) in a complex neighbourhood 
of K 2 . Moreover, the partial derivatives of A with respect to /1 and v are given by 





which yields the following estimates for the growth rate of the eigenvalues (compare [lj, Chap. VII, §3, Sec. 4]): 

< 1. 







dA 








dA 




< 


< 1, 




< 








dv 




dv 



Here, || • || denotes the operator norm of a (2 x 2) matrix. In addition, by [fH Chap. V, § 3, Theorem 4.10], we have 

min |A-A,( K ; 0,0)| < ||T( M , < max{|/i|, \v\} (8) 
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for each eigenvalue A of A(k; fi, v). Finally, by interchanging the components of S(6), we obtain that a point A is an 
eigenvalue of A(n] fi, v) if and only if —A is an eigenvalue of A{— k; /i, —v). Since the eigenvalues depend holomorphically 
on fi and v, the identity 

\j{K;n,v) = -\-j(-n\n,-v) 

holds for all (fi, v) in a neighbourhood of M 2 . Therefore, we restrict our attention to the case k € [i, oo). Note that 
A G C is an eigenvalue of A(k; fi, v) if and only if the system (J5J has a nontrivial solution S(&) satisfying 



I \S(d)\ 2 d9 < oo. 
Jo 



By means of the transformation 



S(0) 



tan • 



cot ■ 



y(sin 2 §), 



e (O.tt), 



the differential equation J5J is equivalent to the system 

y'(x) 



1 

-B 

x 



x-1 



Bi+C 



y{x) 



on the interval (0, 1) with coefficient matrices 



B := 



- 2 ~ 4 A* _ A 







re I 1 
2 t 4 



Bi := 



2 ' 4 

^ — A 



K _ 1 

2 4 



C := 



-2z/ -2fi 
2fi 2u 



and the normalisation condition @ becomes 



Y^— 

I 

X 



y(x) dx < oo. 



(9) 



(10) 



(11) 



(12) 



(13) 



If we consider the differential equation I jllj l for a fixed k 6 (0, oo) in the complex plane, then it has two regular 
singular points, one at x = and one at x — 1 with characteristic values ± (■§ + ¥). From the theory of asymptotic 
expansions (see [24]], for example), it follows that for each A £ C there exists a nontrivial solution 

Vo(x, A) = x% + ^h(x, A), xe^o, (14) 
of I jlljl in the unit disc So C C with centre 0, where h(- , A) : 23q — ► C 2 is a holomorphic function, 



h(x, A) = ^ a: n ft n (A), /i (A) 



fj, — X 



(15) 



n=0 



Here /io(A) is an eigenvector of -Bo for the eigenvalue f + t, and the coefficients h n (X), n > 1, are uniquely determined 
by the recurrence relation 



(Bo-a~ n) h n (X) = (B + B 1 -C+l-a-n) /i„_ x (A) + C/i„_ 2 (A) 



(16) 



with a := f + 7 and /i_i(A) := 0. Since the matrices Bo and £?i depend holomorphically on A, the coefficients 
h n : C — ► C 2 are holomorphic functions. By slightly modifying the proof of [24J, Theorem 5.3], it can be shown 
that the series (j!5j) converges uniformly in every compact subset of QSo x C. Thus, by a theorem of Weierstrass, 
h : So x C — ► C 2 is a holomorphic vector function in the variables (x, A). Now, let 



h (|, A) =: 



/(A) 
5(A) 



and we define the holomorphic function A : C — ► C by 

A(A) :=/(A) 2 - 5 (A) 2 , AeC. 
The following Lemma provides a connection between the eigenvalues of A and the zeros of A. 



(17) 
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Lemma 1 For fixed k £ [i, oo) and (fi, i/) e C 2 , a point X £ C is an eigenvalue of A(k; /i, v) if and only if X is a zero 
of the function A given by (|17l) . This is equivalent to the statement that the differential equation lllfl has a nontrivial 
solution of the form 

y(x) = [x(l-xp+ir,(x), xeC\{0,l}, (18) 

where T] : C — ► C 2 is an entire vector function. As a consequence, if S is an eigenfunction of A(n; /i, v) for some 
eigenvalue X, then 

\S(0)\ < Csm K 8, 0g(O,tt), (19) 

with some constant C > 0. 
PROOF. Defining 

/ n i \ 

(20) 

we have K^ 1 = K and KBqK = B±, KCK = —C. Hence, y is a solution of the system (|11|1 if and only if the function 
Ky(l — x) satisfies <|11[) . In particular, y\(x) :— Ky (l — x) is a solution of (|llf> in the unit disc Q3i C C with centre 
1, and yi has the form 

y 1 (x,\) = (l-x)% + *Kh(l-x,\), xe<Bi. 

Moreover, by the Levinson Theorem (see 0, Theorem 1.3.1]), any solution of (|ll|l which is linearly independent of yo 
in (0, 1) behaves asymptotically like x~^~t. [vq + o(l)] as x — ► 0, where vq is an eigenvector of B$ for the eigenvalue 
— § — j. Similarly, any solution of ljll|l which is linearly independent of y\ in (0,1) has the asymptotic behaviour 

[x — l) - ^ - ^ [v\ + o(l)] as x — * 1 with an eigenvector v\ of B\ for the eigenvalue — 5 — \, Now, if A is an eigenvalue of 
A(k; /i, is), then the system l)llfl has a nontrivial solution ?/ satisfying |JT3J|, and it follows that y(x) = y a (x, X)c a holds 
in (0, 1) with some constants c a 6 C \ {0}, a G {0, 1}. Thus, yo and yi are linearly dependent, and the Wronskian 
W(x,X) := det (yo(x, A), yi(x, A)) vanishes identically for all X G (0,1). In particular, = W {^,X) = 2 _K_ 2A(A). 
Conversely, if A(A) = 0, then W A) = 0, which implies that y$ and y\ are linearly dependent. Hence, yo(x) = y\{x)c 
with some constant c £ C\{0}, and therefore y is a solution of Ijlll) satisfying the condition 113|) on the interval (0, 1). 
Moreover, we immediately obtain that j/o has the form (|18fl with a holomorphic vector function r\ : *Bo U Si — > C 2 , 
and since (|ll|l is regular in C\{0, 1}, we can extend r\ : C — > C 2 to an entire function by the existence and uniqueness 
theorem. Finally, by means of the transformation (|10fl . an eigenfunction S of A(n; [i, v) has to be a constant multiple 
of 



sin - 

sin^( 2 U(sin 2 f), 0€(O,7r), 

cos ■ 



2 / 

and this yields the estimate (|19|l . □ 

Lemma 2 For fixed n £ [5,00) and j £ T, \ {0}, i/ie j-^/i eigenvalue Xj(K;fi,i/) of A(k;[l,v) has a power series 
expansion of the form 

00 

X j ( K ;n,v)= \n, n ti m v n , A 0)0 = A j (k;0,0), (21) 

m,n— 

which is uniformly convergent in the polydisc £ := {(/i, f) G C 2 : |;/| < ^}. Moreover, for all integers m and n, 
the following estimate holds: 

|A m ,„| <(\ K \ + \j\)2 n+m . (22) 

Proof. Since the coefficient matrices in l|llf) depend holomorphically on (A,/i, v) £ C 3 , we can modify |2J, Theorem 
5.3] appropriately in order to obtain that h in (|14fl and therefore A = A(A, [i, v) as given by (|17fl are holomorphic 
functions on C 3 . By a similar reasoning as in the proof of Lemma ^ we can show that for fixed £ C 2 

the eigenvalues of A(k; /i, v) coincide with the zeros of the function A 1 — > A(A, /1, v). In particular for the case 



6 



(/i, v) S R 2 these zeros are simple because A(k; /i, v) has only simple eigenvalues. Hence, by solving the equation 
A(A, fi,v) = and using the implicit function theorem, an eigenvalue Aj(/-c;/i, 7/) of the operator A(k;h,v) depends 
holomorphically on (/i, v) in a complex neighbourhood of R 2 . Furthermore, the estimate (jHJ implies that the set 
{A € C : miiij^o |A — Xj(n; 0, 0)| > 5} contains no eigenvalues of A(k;h,v) for all v) G £. Thus there exists a 
holomorphic solution A : £ — > C of the equation A(A, /i, v) — 0, which is uniquely determined by A(0, 0) = Aj (k; 0, 0). 
Consequently, Aj(«; /i, z/) is holomorphic in £, and therefore it has a power series expansion in £ of the form (|21|l . In 
addition, by Cauchy's formula, 

Xj(K;iJ,,u) 

71 + 1 ~ ' 



47T 2 J ae /Li m+1 J/ 



and applying (|SJ) and JTJl, it follows that 



\Xj{K)^v)\ < |A i (K;0,0)|+max{|/x|,|^|} < \k\ + \j\ 



□ 



which gives the estimate (I22|) . 

According to Lemma ^ f° r fixed parameters v) E C 2 the eigenvalues of A(k; /i, v) are exactly the zeros of the 
function A(A) given by l|17|l . In principle, this result can be used for numerical computation of the eigenvalues. 
However, in order to calculate A(A) at some point A £ C, we first have to determine the coefficients h n (X) with the 
help of the recurrence relation and subsequently we need to evaluate h(x, A) at x = ^ by means of the power 
series expansion (|15|) . Unfortunately, this method requires the calculation of two consecutive limits, making things 
rather complicated. In the remaining part of this section we show that there is yet another function O which encodes 
the eigenvalues of A(k] h, v). The main advantage of O is, that it can be obtained by only one limit process. 

By setting y(x) :— x a (l — x) 1 ~ a y(x) with a := -| + 4, the system ljTT)l becomes 



y'{x) = 



-S + — 

X x - 



1 



B l +C 



with the coefficient matrices 



Bn := 



— K —5 [1 — A 



, B 1 := 



-i 



°/ v^ _A_1 

Now, there exists a holomorphic solution of ill' Ml) in Q5i given by 



, c = 



(23) 




y(x,X) = J2x n d n (X), d (A):= 



71=0 



(24) 



where do(X) is an eigenvector of Bq for the eigenvalue 0. In addition, the coefficients d n (X), n > 1, are uniquely 
determined by the recurrence relation 



d n (X) = {B - n)- 1 [{E - n)d n -x{X) +Cd n - 2 (X)] 



with 




Finally, we denote by B„(A) the second component of d n (X). 

Lemma 3 Let n € [i, 00) and (/i, v) € C 2 be fixed. Then, for each X G C, </ie /«mi£ 

6(A) := lim 8„(A) 



(25) 



exists, and Q : C > C is a holomorphic function. Moreover, a point X G C is an eigenvalue of A(k; fi, v) if and only 

'/(-)■: A) 0. 
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Proof. For fixed A € C, the differential equation (12,'j[| has a regular singular point at x = 1 with characteristic values 
— f and k — \. First, let us assume that their difference k + j is not an integer. In this case the system (12. '{f has a 
fundamental system of solutions in a complex neighbourhood of x — 1, which can be written as 

oo oc 

yx(x, A) = (1 — x)- 1 £(1 - x)"4(A), y 2 (x, A) = (1 — £(1 - x) n d£(A), (26) 

?;.=0 ?;.=0 



where 



0\ „ ( K+h 



4(A) = I =: e 2) 4(X) 



2 



x 1 / \^~ A / 

are eigenvectors of B\ for the eigenvalues —1 and k — |, respectively. Now, j/ can be written as a linear combination 

y{x,X) = ji(X)yi(x,X) + 72(A)y 2 (£, A) 

with connection coefficients 71(A), 72(A) £ C. Applying [2^, Corollary 1.6] to the system (12.'S[) gives 

lim d n (X) = 7 i(A)e 2 , (27) 
n — >oo 

and therefore the limit (|25|) exists. Furthermore, A is an eigenvalue of A(k; fj,, v) if and only if 71(A) = 0, i.e., if and 
only if 0(A) becomes zero. Finally, it can be shown that the functions d n converge uniformly in every compact subset 
of C, and Weierstrass' theorem implies that is an entire function. 

Now, suppose that k := k + ^ is a positive integer. In this case, a fundamental system of the form l|26|l may not 
exist. Nevertheless, it can be proved (see Lcmma[f)]in Section VI) that the system (|23f) has a fundamental matrix 

00 

Y(x, A) = G(A) ]T H n (X)(l - x) n {\ - x) D {l - x) J ^\ 

n=0 

in a complex neighbourhood of x = 1, where D := diag (—1, k — 1), Hq(X) — I and 

G(A) H^-a)' J(A) =U») 

with some q(X) G C. In particular, we can write Y in the form 

00 

Y{x, A) = H(x, A)(l - x) J M, H{x, A) = £(1 - x) n D n (X), 

n=0 

where 

/ o o \ - / -1 \ 

Since y solves the system (|23H . there exists a vector c(A) S C 2 such that y(x, A) = F(x, A)c(A), and [2ll Theorem 1.1] 
implies 

d„(A) = D (A)^ (-J(A)) r(n + 1)^ (n - J(A)) c(A) + O (n 5 - 1 ) (28) 

for arbitrary <5 > 0. For the definition and discussion of the reciprocal Gamma function for matrices we refer to the 
Appendix in plj . Particularly, for the Jordan type matrices — J(A) and n — J(A) i we obtain 

1 / =„.\ / 1 0\ 1 / ~,.,\ / r^T) \ 



(-J(A)) = f : , ±(n-i(A)) 



rV 7 V* V ' r V 7 V * rRI) / ' 

Now, if 71(A) denotes the first component of c(A), then (|28() implies l|27ll . Since A is an eigenvalue of A(n; /i, v) if and 
only if 71(A) = 0, the proof of Lemma|21is complete. □ 
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III. A PARTIAL DIFFERENTIAL EQUATION FOR THE EIGENVALUES 

Theorem 1 For fixed k£ |, oo) and j € Z \ {0}, the j-th eigenvalue A = Xj(n; fi, v) of A is an analytical function 
in (p, v) G M 2 satisfying the first order quasi-linear partial differential equation 

dX dX 

(u-2v\) — + (v -2uX) — + 2ku+2ui> = 0, (29) 
o/j, ov 

where Aj(K;0,0) is given by (|T|). 
Proof. Let 

/ Si(0) \ 

m = : U<») ) ■ ^ ' e(0 •' , • 

be that eigenfunction of A(k; /x, v) for the eigenvalue A = Xj(n; fi, v) which is normalised by the condition (S*, S) = 1. 
Introducing the functions 

U{9) :=S 1 (6) 2 + S 2 (9) 2 , V(9):=S 2 (9) 2 -S 1 (9) 2 , W{9) := 2 Si(0) S 2 (9), 

a straightforward calculation shows that U, V, and W are solutions of the system of differential equations 

U'{9) = 2 ( v sin 9 + -?—) V(9) + 2ficos9 W(9), (30) 

V sin 9 1 

V'(0) = 2(v sin 9 + -^—) U(9) + 2 X W{9), (31) 

V sin 9 J 

W'{9) = 2[icos9U(9) -2XV{9). (32) 
Now, from analytic perturbation theory (compare [lj, Chap. VII, §3, Sec. 4]) it follows that 
dX fdA 



■ )n \ )fl S ' S l = l S ( 9 Y \ n n )S(9)dO=l cos9V(9)d9, (33) 



cos 9 








cos 9 





— sin 


sin 9 






^.= (^S,S)= I S(9)*\ , , ; )S(9)d9 = -l sin9W(9)d9. (34) 



dv \ dv 

In addition, from 1|19|) we obtain the estimates 

\U(9)\, |V(0)|, \W{9)\ <Csin 2K 9 

with some constant C > 0. Since k is positive, U, V and W vanish at 9 = and 6* = tt. If we integrate i|33|) by parts 
and replace V'(9) with the r.h.s. of l|31|l . then we get 



dX 

Of l 



Taking into account that 



we have 



/ sm9V'(9)d9 = - (2vsm 2 9 + 2 K )U{9) + 2Xsm9W(9)d9 
Jo Jo 

-(2v + 2k) / U(9)d9-2X / sinfl W{9) d9 + 2 v / cos 2 9 U{9) d9. 



/ U(9)d9 = (S,S) = l, / sm9W(9)d9 = -— , 
Jo Jo ov 



dX dX 
H— = -/j (2v + 2 k) + 2ijX— +2/JV / cos 2 9 U(9) d9. 
dfi dv J 

Moreover, equation 1(3*2^) implies 

2 (i cos 2 9 U{9) = cos 9 W'{9) + 2 A cos 9 V(9), 



(35) 
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and integration by parts gives 

[>TT I>7Z />7T 

2y,v I cos 2 6 U{9) AO = v I cos<9 W'(6) d9 + 2 v A / cos 9 V (9) d9 
Jo Jo Jo 

= v \ sin9W(9)d9 + 2vX cos 9 V(9) dd = -v — + 2 v A — . 
Jo Jo dv 



(36) 



Replacing the last term on the r.h.s. of (|35|l with (|36|1 . we obtain exactly the partial differential equation (|29|) . □ 

The PDE (|29|l can be used in order to derive a power series expansion for Xj with respect to fi and v. For this 
purpose we introduce the new coordinates (compare psT ]) 

a := v — fi, f3 := v + /i- 
Then X(a, (3) := Xj ^k; ^§^J is a solution of the transformed partial differential equation 

a(l + 2A)^+/3(l-2A)g = K (a-/3) + i( a 2 -/? 2 ), (37) 

where A(0,0) = Aj(re;0, 0) is given by Q- As A depends analytically on (a, (3), there exists a series expansion for A 
of the form 

oo 

X(a,(3) = cm, n a m p n (38) 

m,n— 

(for clarity, the indices re and j in the coefficients c mjIl and in the function A have been omitted). Furthermore, 1371) 
is equivalent to 

fax ax 2 \ n (dx dx 2 \ , „. i, 2 n2 . 



da da } \ d(3 8(3 



and since 



oo / m n 



A(a,/3) 2 = £ p^c r , s c w aT 

m,n=0 \r=0 5—0 / 



we obtain the identity 



^ (m + rc)c m ,„ + (to - n) ^ ^ c r , s c m _ r ,„_ s ) ot m (3 n = k(o - (3) + - (a 2 - /3 2 ) . 

m.n—0 \ r—0 s—0 / 

Comparing the terms of equal order in a and (3, it follows that 

co.o = Aj(k;0, 0) =: c , c 10 = - —— , c .i = - ; 

2 Co + 1 2 Co — 1 

_ (2c + l) 2 -4 K 2 _ n _ (2c -l) 2 -4 K 2 

C2 '°~ 4(2 co + 1)3 ' Cl - 1 " ' C °' 2 " 4(2 c -1)3 ' 

and for to + n > 2 the coefficients c m .„ satisfy 

((m + n) + 2co(m- rj))c m ,„ = (n- m) c r!S c m _ ri „_ s (40) 

(r,s)e[m,n] 

where [to, n] denotes the set of all pairs (r, s) £ Z 2 with < r < to, < s < n and 0<r + s<m + n. In particular, 
if to = ri > 0, then we get 2 n c n/n = 0, which implies 

c, v „ = for all n > 0. (41) 
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Moreover, if re is not a rational number, i.e., re G [i, oo) \ Q, then the initial value Co is not a rational number, and we 
have (to + n) + 2 co(to — n) ^ for all (m, n) € Z 2 with m + n > 2. In this case (|40|) gives a recurrence formula for 
all coefficients of the power series expansion (|38ll . 

Now, suppose that k£ [i oo) n Q. Then c is a rational number with |c | > 1, and we get f^fpj = | with some 
coprime integers p and g. Now, the prefactor on the l.h.s. of 140f> becomes zero if and only if to — £p, n — Iq with some 
positive integer £, and thus the coefficients ci p ^ q are not determined by l|4()(l . However, we can by-pass this problem 
if we regard re as an additional parameter in our eigenvalue problem. Since the coefficient matrix of f|l 1 1) depends 
holomorphically on re G C + := {z G C : Rez > 0} and (A, p, v) G C 3 , we obtain in a similar way as described in 
Section 1 that h in (|14f> and therefore A = A(re; A; //, v) given by (|17|l is a holomorphic function on C + x C 3 . Moreover, 
in the same way as in the proof of Lemma^] we can show that for fixed re G oo) and (/x, v) G C 2 the eigenvalues of 
A(k; /i, v) coincide with the zeros of the function A i — ► A(re; A; /i, v). In the case (re, fi, v) € [h, oo) x K 2 these zeros are 
simple, since A(k; /i, v) has only simple eigenvalues. Hence, by solving the equation A(re; A; /x, v) — 0, we find that an 
eigenvalue Aj(re; /i, v) is a holomorphic function in a complex neighbourhood of [A, oo) xl 2 . In particular, A depends 
holomorphically on (re; a, /3), and for a given re € [5, 00), there exists a power series expansion of the form 

00 

A(re + e; *,/?)= ]T c«„ £ <a m /T 

/,m,n— 

in a neighbourhood of (k, 0, 0). In the following we derive a recurrence relation for the coefficients Cm,n- Since 

oo / I m n \ 



A(re + E ;a,/?) 2 = £ (EEE^W')^ 

from l|39|l it follows that 



.(*) ,.(«-*) 

l.m,n=0 \t=0 r=0 s=0 



/ m n 



£ (m + n)c« n + (m - n) £ £ £ c« c^_ £ <c/"/?» = re(a - /?) + e (a - /J) + ~ (a 2 - /3 2 ) . (42) 

l,ra,n=0 \ t=0r=0s=0 / 



Moreover, Q implies that 



sgn(i)(re-i + |i|), if 1 = 0, 
c o.o = 7TT r7 («;0,0)= i sgn(j), if Z = 1, 



II 6W 



0, if Z > 1. 



Comparing the terms of equal order in 14211 . we obtain 



c (°> - A fre- 01 - m c (0) - K c (0) - K 

c 0,0 — A J \ K > u ' u i — • c 0i c l,0 — o 11' 0.1 — o 1 ' 

Z Co + 1 Z Co — 1 

(0) _ (2c + l) 2 -4re 2 jo) = Q jo) _ (2 cq- l) 2 -4re 2 



2 <°~ 4(2c + l) 3 ' M ' °' 2 ~ 4(2 c -l) 3 
(1) _ 2 c + 1 - 2sgn(j)re (1 ) _ 2 c - 1 - 2 sgn(j) re 
Cl >°- (2c + l) 2 ' C °' 1_ (2c -l) 2 ' 

while the remaining coefficients are determined by the identity 

((m + n) + 2co(m-n))cW n + (m-n) £ ^0^ = 0, Z + m + n>2. (43) 

(i,r,s)e[i,m,n] 

Here [Z, to, n] denotes the set of all triples (t, r, s) G Z 3 with 0<i<Z, 0<r<m,0<s<n, and < t + r + s < 

I + to + U. In the case (to + n) + 2 co(to — n) — 0, the prefactor of Cm,n in (|43|l vanishes, and since m — n =^ 0, we get 
for Z > 

0= V c^c^ -2c (1, c (, - 1) + V c (t) c ( ' _t) (44) 

(l,r,s)e[!,m,nl (t,r,s)ep,ro,n]* 
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where [Z, to, raj* := [Z, m, n] \ {(1, 0, 0), (Z — 1, m, n)}. Now, for all coefficients Cm, n with Z + to + n > 2, ((43(1 implies 



(•;;;„ = ; ^ 4*s c m-r,n-s if (m + n) + 2ca(m- n) ^ 0, 



(m 

whereas (l4*4*)l and c^q = sgn(j) yield 



fm + n ) + 2 cn(m — n) 

(t,r,s)e[l,m,nl 



= _ sgnO) ^ c (,) c a-^_ s jf ( m + n ) + 2co (m - n) = and Z > 1. 



2 

(t,r,s)e[Z,m,ra]* 



These recurrence relations can be used to determine all the coefficients c rn . n — Cm, n of the power series expansion i|38|) 
in the case that k is a rational number. 

Remark 1 A series expansion for the eigenvalues A with respect to (a, (3) has been given by Suffern, Fackerell & 
Cosgrove 123. Sec. 8], however, only the coefficients c m . n with m + n < 5 have been determined. Furthermore, Kalnins 
& Miller flSy studied a series expansion A — X)^Lo ^« a ™ f or the eigenvalues in terms of the Kerr parameter a, but 
also in this paper only a finite number of coefficients Ao, . . . , A3 have been explicitly computed. A general recurrence 
relation for the coefficients of (I38|) could not be found in the literature. Moreover, the problem of dividing by numbers 
which may be zero has not been noticed in \2oV and \1S\L Finally, it should be noted that some of the diagonal entries 
c n ,n for n > in '28, Table I] are not equal to zero, in contrast to our result 141(1. 

IV. SOLUTION OF THE PDE BY THE METHOD OF CHARACTERISTICS 

In this section we study the PDE (|29|l for real parameters (jx,v) € R 2 and fixed k € [5,00) by the method of 
characteristics. In particular, we obtain an exact formula for the eigenvalues in the case = \u\, and for \fi\ =/= \v\, 
it turns out that the characteristic equations can be reduced to the third Painleve equation. 

Theorem 2 Let k € [h, 00), j E 1 \ {0} and r e {— 1, +1} be fixed. Then 



\j(n;n,Tn) = ~ +sgn(j)^ (a.,-(k; 0,0) - +2tk/j + /i 2 , 



(45) 



where 0, 0) = sgn(j) (k — | + |j|) . In particular, if j = t, then 

Aj(«;;/i,T/i) =r f K + ^J + l L - 

Proof. According to Theorem the function X(fi,u) := \j(n\n,v) solves the partial differential equation 1(29(1 . 
Defining u>(/i) := A (/i, r/i), fj, G R, for some fixed r S { — 1, +1}, we obtain 

«j'(aO = 0"> T v) + ^' ' 
and with the help of (|29|) it can be shown that 

liw'(iT) — 2 r/z w(n)w' (/i) = —2 K/i — 2t/U 2 . 
Dividing the above differential equation by —r/i and integrating gives 

- = C + 2rKit + /i 2 , /16I, (46) 
where the constant of integration C is uniquely determined by 

0= (w(0)-^Y = (A,( K ;O,O)-0 2 . 

Now, from ((46(1 it follows that 

wQjl) = T -+e ^(A J ( K ;O,O)-0 2 + 2r KAt + ii 2 (47) 
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with some e £ { — 1, +1} and the square root assumed to be non-negative. We have to take the sign e such that the 
l.h.s. of (|47|) is analytic and coincides with Xj(k; 0, 0) at the point /i = 0. If j = r, then Xj(k; 0,0) = r («+ 5), and 
l|47Jl implies «;(//) = J + £ (tk + /j). Inserting /1 = 0, it follows that e = r, i.e., = r (ft + i) + /i. In the case 

j 7^ r we have I (re; 0, 0) — J I > re + 1 and thus the radicand in (|4T|> is positive for all \i £ R. Moreover, by means of 



A i (rt;0,0) = ^(0) = \+£\J (Aj(k; 0, 0) - I)* = I + e | A,-(k; 0, 0) - 
and (J7J, we get e = sgn (A.,(re; 0, 0) — ?) = sgn(j), which completes the proof. □ 



Remark 2 For a given half-integer re and (i — v, this result has been shown by Suffern, Fackerell & Cosgrove using a 
power series expansion for the eigenfunctions of QJ - (J2J in terms of hypergeometric functions (see J23 . Sec. 3-5]). 
Here, we obtained the formula for Aj(re; fi, ±/i) as an immediate consequence of the partial differential equation l|29|) . 
Moreover, it should be noted that the formula 0, (54)] given by Chakrabarti is not correct. 

Now, let us consider the case |/x| ^ \v\. For this reason, we introduce new coordinates (t, v) £ (0, 00) x (M \ {0}) by 

fl (t,v)= t -(y + ^j, v {t,v)= t -(v-^) (48) 

with some fixed a £ {— 1,+1} (note that a — ±1 corresponds to the cases \n\ > \v\ and \n\ < \u\, respectively; 
moreover, this transformation maps v — const, onto lines in the (/i, ^)-plane starting at the origin, while the curves 
t = const, are mapped onto hyperboles). By setting w(t, v) = A(/i, v), we have 

dw If o-\d\ If <r\ dX 1 / dX dX 



and H29|l becomes 



dt 2 V" ' vJ dfi ' 2 V v S dv t \^ dfx dv 
dw t f ' a \ dX t f ' a \ dX 1 / dX dX 



dv 2 V v 2 J dfi ' 2 V v 2 J dv v \ dpi ^ dv 



dw 2vw dw f a \ t ( 1 



dt t ft, + H tt+ y + iK-^)= - (49) 



The characteristic equations of this PDE are given by 



^ ) = -«(„( t ) + _£ y )-^(„ (t ) 2 -_i p ). ,51, 



From (|50|l we obtain that w(t) = — fU , and l|51|) implies 



2v(t) 

v'(t) tv"(t) tv'(t) 2 _ ( a \ t ( 2 1 



"(*) + 77^ +o " 



2v(t) 2v(t) 2v{t) 2 V v(t)J 2 V «W 2 

Multiplying the above differential equation with 2v(t) 2 , we get the following third Painleve equation 

tv v" - t(v') 2 + v v' - 2k (v 2 + a) v - t (v 4 - l) = 0, (52) 

with parameters a = a[3 = 2k and 7 = —S — 1 (see, for example, 01 or For further details on the Painleve III 

we refer to e.g. 0, HJ and @. 

In general, Painleve III is not solvable in terms of elementary functions, and therefore we cannot expect a closed 
expression for the eigenvalues of A(k; fi, v) in the case |/x| 7^ \v\. On the other hand, for particular values of k there 
exist so-called special integrals of polynomial type for this equation, i.e., polynomials Q in t, v and v' with the 
property that every solution of the differential equation Q(t,v,v') = satisfies As it will be shown below, such 

special integrals are related to algebraic solutions of the PDE (|29|) . i.e., solutions, which are zeros of a polynomial in 
A with rational coefficients in /i and v. Moreover, taking into account that the eigenvalues Xj(K; fi, t/j,), t £ {—1, 1}, 
of A(k; fi, Tfi) satisfy the quadratic equation 

(x-^ 2 = C + 2tkii + with C := (Xj{k;0,0) - , 

the question arises if such an algebraic expression for the eigenvalues of A(k; fi, v) exists in the case \n\ 7^ \v\. A first 
step towards the answer of this problem is given by the next Lemma. 
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Lemma 4 Suppose that there exists a polynomial 

N 



P(A;/i^) = J]P„(M^)A n , P N = 1, 



n=0 

of degree N > in A with rational coefficients P n in fi and v such that the zeros Zj(jx, v), j = 1, ... N , of P{ ■ ; fi, i>) are 
simple, and that the functions A = Zj are solutions of the partial differential equation H29|) . Then k is a half-integer. 
Moreover, if N — 1 or N — 2, then k = \ and P(X; fi, v) — (A + ft) . 

PROOF. Let Q(t, v, v') — Yl n =o Qn(t, v)(v') n be the polynomial in v' with coefficients 

2 v 



,(«,«):= I-—) P n (ji(t,v),v{t,v)), n = 0,...,N. 



where fi(t, v) and v(t, v) are given by gSJ. Note that the Q n are again rational functions in t and v. Moreover, let 
v : D — ► M\ {0} be any solution of the first order ODE Q(t, v, v') = on some interval £> C R\ {0}. For the function 

te ®> (53) 

we obtain 

= Q (t, v(t),v'(t)) = J! Qn (t, v{t)) (- 2 ^ } t W(t) ^ n = (-^) N P («»(*);/*(*, "(*)), «(*))) , 

and thus is a zero of P (■ ; fi (t,v(t)) , v (t,v(t))) for each t € 35. Since the zeros of this polynomial depend 
analytically on the parameter t according to the implicit function theorem, there exists an index j € {1, ... , N} 
such that w(t) = Zj (fi(t,v(t)) ,v (t,v(t))) for all t e S. Furthermore, as Zj solves the PDE (|29|l . it follows that 
(t,v(t),w(t)), t E D, is a characteristic curve of Q49J1 . and thus w is a solution of l|52|l . Hence, Q(t,v,v') = 
implies (|52|) . and therefore Q is a special integral of rational type for this Painlcve III. Multiplying Q(t,v,v') by 
an appropriate polynomial r(t,v) in t and v, we obtain that Q(t,v,v') := r(t,v)Q(t,v,v') is a special integral of 
polynomial type of degree N with respect to v'. Now, Theorem 2] yields that such a special integral exists if and 
only if 2 k ± 2 ctk = 2(2/c — 1) with some integer k, i.e., k = k — i is a half-integer. In addition, by 16, Lemma 3], the 
relation (aq — p)n = N has to be satisfied for some integers p, q G {— N, —N + 2, . . . , N — 2, TV}. In the case N = 1 
or N = 2, these conditions imply k = ^, and the corresponding special integrals of polynomial type are explicitly 

known, namely r(t) v s (y' + v 2 + o~) , where r is some polynomial in t, and s is an integer (compare |lfil Section 2]). 
Hence, Q(t, v, v') = (v' + v 2 + cr^j N and 

/ t \ N 

P (w(t);n(t, «(t)), u(t))) = ( - w^J ) (*. W W) = ( W W + ^ ((*. «(*)))) - i G 2), 

which yields P(A; /Lt, t/) = (A + m) W if N = 1 or N = 2. □ 

As a consequence of this Lemma, if a solution X(jjL,v) of the PDE 129|) is a zero of a linear or quadratic polynomial 
with rational coefficients in /z and v, then k = i and A(/x, j/) = — /x. In fact, the function A(/x, i/) = —fj, solves (|29H for 
/t = i, but since A(0, 0) = and the spectrum of A(h; 0, 0) is given by Z \ {0}, it is not an eigenvalue of A(h; /x, v) for 
any (/i, v) el 2 . The following considerations show that this solution is nevertheless of interest. 

V. MONODROMY EIGENVALUES 

In this section we consider the case that k is a positive half- integer, i.e., k — k — ^ with some positive integer 
k, and we assume that the matrix C defined in (|12|l has distinct eigenvalues, i.e., /x 2 7^ v 2 . For such k and (/x, v) 
there is in addition to the classical eigenvalues of A(k; /x, v) another type of "special values" which we call monodromy 
eigenvalues. In order to introduce this concept, we first recall the characterisation of eigenvalues according to Lemma|] 
A point A is an eigenvalue of A{k\ fi, v) if and only if the system {Tj has a nontrivial solution of the form 

y(x) = [x(l - x)] ! i + ^ri(x) (54) 
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where r\ : C — ► C 2 is an entire vector function. Now, as the difference of the characteristic values ± (§ 
1 is an integer, the differential equation Hll|) has a fundamental matrix of the form 



at and 



Y(x) = [x(l-a:)]" 



(55) 



where H(x) = H a (x)(x — a) Ja in 5B Q , a 6 {0, 1}, with some holomorphic function H a ■ 93 a — ► M2(C) and a Jordan 
matrix J a (see |24l Theorem 5.6]). Hence, the matrix function H is in general not holomorphic in 23o U 2$i since it 
involves logarithmic terms. In the following, a point A S C is called monodromy eigenvalue of ^4(k; //, v) if and only 
if the system l)llfl has a fundamental matrix of the form (|55|l with the property that H : C — ► Ma(C) is an entire 
matrix function. Monodromy eigenvalues are characterised by the following Lemma. 

Lemma 5 For a given half-integer k > and (/z, v) G C 2 with /i 2 ^ v 2 . a "point A G C is a monodromy eigenvalue of 
A(k; /i, v) if and only if the system l|ll|) has a nontrivial solution of the form 



[x(l - x)] ^ 3p(x)e 2 



(56) 



where p : C — ► C 2 is a polynomial vector function and t = ±\J v 2 — . 

Proof. By means of the transformation y(x) = x a (l — x) a y{x) with a 
equivalent to the system 



the differential equation l|ll|) is 



S'lx) 



1b 

X 



1 



(57) 



where 



Bn 




Bi 




(58) 



Now, if A is a monodromy eigenvalue of A(k;/j,, v), then the system (|57|) has a holomorphic fundamental matrix 
H : C — ► M 2 (C). Since the coefficient matrix of l|57|l is a rational matrix function which is bounded at infinity, 
an extension of Halphen's Theorem (see 0, Theorem 2.4]) implies that the system (|57f) has a fundamental matrix of 
the form R{x)e Dx with some rational matrix function R and D := diag (— 2t, 2t) (note that ±2t are the eigenvalues 
of C). Moreover, R(x)e Dx — H(x)Q with some invertible matrix Q, and therefore R(x) — H(x)Qe~ Dx is an entire 
matrix function in C. This implies that R : C — ► M 2 (C) is a polynomial. Vice versa, suppose that the system 111|) 
has a nontrivial solution y(x) = [x(l — x)]~ a p{x)e 2tx with some polynomial vector function p : C — ► C 2 . Defining 



-21 



Ky(l -x) = [x(l - x)]- a Kp{l - x)e 



-2tx 



with K given by (|20|l . it follows that y is a solution of Ullfl which is linearly independent of y. Therefore, Ijlll) has a 
fundamental matrix of the type (|55|l . where H(x) — (p(x)e 2tx , Kp(l — x)e~ 2tx ) is an entire matrix function. □ 

Theorem 3 For fixed k — k — ^ with a positive integer k there exists a polynomial P(k; A; fi, v) of degree 2k — 1 in 
A with polynomial coefficients in /i and v such that for each (/i, v) £ C 2 with ^ v 2 a point A G C is a monodromy 
eigenvalue of A(k; /z, v) if and only if X is a zero of P{k; ■ ; fi, v). Moreover, the integers 1 — k, . . . , k — 1 are the zeros 
of P(k; ■ ; 0, 0), and for k — | we obtain P(h; \;n,v) = A + /i. 

Proof. A point A is a monodromy eigenvalue of A(k; /i, v) if and only if the differential equation l|57|l has a nontrivial 
solution p(x)e 2tx , where p(x) — X)^=o^" a; "' Pjv ^ 0, is a polynomial vector function, and t — ±^v 2 — /i 2 . In the 
following we assume t — y/ v 2 — (the main branch of the square root) but all considerations remain valid if we 
replace t with —t. If we set A := A — /i and C := C — t, then the polynomial p satisfies the differential equation 



p'(x) = 



x 



Bi+C 



(59) 



where the coefficient matrices take the form 

'0 -A\ 

k ' 



Bn = 



B 1 = 




C = 



-2v-2t -2/x 
2fi 2v - 2t 
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It is easy to see that the coefficients p n S C 2 , n — 0, . . . , N, form a nontrivial solution of the linear system of equations 

B oPo = 0, (Bo - l) P i + S Po = 0, (60) 

(B -n)pn + {S + n-l)p n -i-C Pn -2 = (n = 2,...,N), (61) 

(S + N) PN -C PN ^ = 0, -C PN = Q, (62) 

where 



-2v-2t-k -2/i + A 
2^i + A 2v-2t-k 



S : — C — Bq — B\ — 

Multiplying the first equation in (|62|l from the left with the matrix C + 4t and observing that (C + 4t)C — 0, we get 

{ N-k -A \ „ 

= (C + At)(S + N)p N = [ )Cp N + 4t(N -k)p N = 4t(N ~k)p N . 




Since t ^ and pm ^ 0, it follows that N = k. Due to technical reasons we have to distinguish between the cases 

k > 2 and k = 1. We will proceed at first with a detailed proof for the more complicated case k > 2. Adding the 
second equation in (|62|l to the first one and then both equations in (|62|l to H61(l for n = N, we obtain 

B oPQ = 0, (B - l) Pl + Sp Q = 0, (63) 

(B o -n)p„ + (S + n-lK-i-Cp„_ 2 = (n = 2,...,k-l), (64) 
/ -1 A \ 

-Bipk + ( J p*-i - Cp fc _2 = 0, (65) 




Pk - CVk-i = 0, -Cp k = 0. (66) 

The system above consists of 2k + 6 linear equations for 2k + 2 unknowns. In the following we prove that only 2fc + 2 of 
these equations are linearly independent. Summation of all equations l|6(J [) - l)62|) yields —B\ Y^=oP n = 0- Because of 
rank (B±) — 1, it is possible to eliminate the second line of the first equation in (|65() by means of line transformations, 
and since also rank (So) = 1, we can delete the first line of the first equation in (|63|l . Thus, the system (|63|) l|66f) 
consists of at most 2k + 4 linearly independent equations. In order to reduce the equations (|66|) further, we have to 
consider the cases v — t ^ and v — t = separately. First, let us assume that v — t ^ 0. Multiplying the equations 
in H66|) from the left by the invertible matrix 



T 




it follows that H66[) is equivalent to 





P*-i=0, I „ „, \Pk=0- (67) 
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Now, we can represent the system of the linear equations (|63|l - (|65|l . (|67|l as a matrix equation Tp 
p= (po,---,Pk)£ € 2k+2 and the (2k + 2) x (2fc + 2) matrix 



with 



/ 



r := 



k 



















q 


-DO 


J 

— 1 

















u 


—C 


s 


f z 


So 


-2/ 













* 







c 


5 


f 2/ 


B -31 

























-C 


S^ 


- (k - 2)1 


B -(k 


-1)1 



















2v 


+ 2t 2(i 


-1 


A 


-k 





















-2/i 2t 


- 2;/ 


A 



























-2/i 2t-2v 



\ 



(68) 



Let f be the (2k + 1) x (2k + l)-matrix obtained from f by deleting the last row and column. Then A is a monodromy 
eigenvalue of A(k; /z, v) if and only if = det T = (2t — 2v) det T, i.e., det T = since v — t ^ 0. Now, suppose that 
v — t = 0. We will prove that also in this case A is a monodromy eigenvalue if and only if detf = 0. Note that v = t 
implies /j = 0, and therefore the equations in l|66(l are equivalent to 




Pk-i = 0, 




(69) 



If A is a monodromy eigenvalue of A(n; fi, v), then the vector p — (po, ■ ■ ■ ,Pk) is a nontrivial solution of the matrix 
equation Tp — even though T is not invertible for v — t = 0. If we assume detT ^ 0, it follows that A ^ and 
the first 2k + 1 components of p must be zero. In particular, pi = ■ ■ ■ = Pk-i = 0, and the first equation in 1691) 
yields pk = 0. Thus p = 0, and this contradiction implies detT = 0. Conversely, if detf = 0, then either A = and 
p := e2k+2 (the 2k + 2-nd unit vector in C 2k+2 ) is a nontrivial solution of (|63|) - (|65|l and l|69|l . or A ^ 0. In the latter 
case, there exists a vector p ^ with components pi, . . . ,f>2k+i & C such that Tp = 0. Defining q :— jrft2k-l, then 
p := (p, q) 6 £?k+2 j g a non t r ivial solution of the equations (|63|l - l|65|) and (|69|l . i.e., A is a monodromy eigenvalue. 
Hence, we have shown that for all (/z, v) £ C 2 with /i 2 ^ v 2 a point A £ C is a monodromy eigenvalue of A(k; fx, v) 
if and only if detT = 0. In order to prove that detT is a polynomial in A of degree 2k — 1, we apply once more 
appropriate line transformations to T. Adding successively the second to the fourth line, the fourth to the sixth line 
and so on up to the 2fc-th line, then detT = detT with the (2k + 1) x (2k + l)-matrix 



k 




















So 


B -I 











—R 


Si 


Bo -21 








Q 


-R 


5 2 B -31 









Q 


Q 


—R 5*3 


B -41 






Q 




Q 


—R Sk-2 B - 


-(fc-l)J 





Q 






Q Q 


—R 


-k 






A 



where 




R := 




-2v -2t-k -2/i 
2/i + A 2v - 2t + n 



,fc- 2. 
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Now, A appears at most once in each row and each column, whereas only the first and the 2k-th line contain no 
entry involving A. It is easy to verify that det T(A; fi, v\ i) has the form ±A; 2 A 2fe_1 + (terms of lower order in A), 
and therefore det T(A; n, v; t) is a polynomial in A with polynomial coefficients in /i, v and t. Moreover, for all 
(/i, v) € C 2 with /j 2 ^ v 2 a point A is a monodromy eigenvalue of A(k; /i, v) if and only if the determinant of 
A; /i, v; t) vanishes. As mentioned at the beginning of the proof, this result remains valid if we replace t with 
—t. Hence, the zeros of the polynomials det A; /i, v; t) and det T(k; A; /i, v; — t) coincide, which implies that 
det T(k; A; fi, v; t) = det A; /i, v; — t). Consequently, the polynomial P(n; A; /i, v) := det T(k; A; /i, v; t) contains 
no terms in t of odd order, and the terms of even order in t depend only on t 2 = v 2 — ii 2 . It follows that P is a 
polynomial of degree 2k — 1 in A with polynomial coefficients in \x and is, and the zeros of P are exactly the monodromy 
eigenvalues of A(n; js, v). 

Next, we prove that the integers 1 — k, . . . , k — 1 are the zeros of the polynomial P(k; A; 0, 0). To this aim, let T 
be the (2k x 2k) matrix obtained from T for (/i, v) = (0, 0) by deleting the last row and column. Then 



Ok 







\ 



Q B -I 







Q + I B -21 
Q + 2I B -2.I 





V o o 



Q + (k- 3)1 B -(k- 2)1 

Q+(k-2)I B -(k-l)I 

-1 A 




and detf = A det IV Moreover, det To = if and only if the equation ^o(p n )n=o = nas a nontrivial solution. Such 
a nontrivial solution is a constant multiple of the vector given by the recurrence formula 



Po 




p n = (n-B ) 1 (Q + n- l)p„_i 



1 



n(k — n) \ 

for n = 1, . . . , k — 1. By induction, it can be shown that 



(k - n)(k + 1 - n) - X 2 
nX 



n-l 



Pn=(-iT [k 1 1)! n[( fc -j) 2 -A 2 ] 



J'=l 




Pn-1 



Multiplying the vector (p n )n=o fr° m the left by the last line of To, we get 

k-2 



0= -] 



("I)*" 1 Jkhy, II - - A2 ] i- k + A 2 + (fc - DA 2 ) = n [(* - jf A 2 ] . 

Hence, det To = if and only if A 2 £ {1, . . . , (k — l) 2 }, and therefore 1 — k, . . . , k — 1 are the zeros of P(k; ■ ; 0, 0). 

It remains to deal with the case k = 1, where we have to consider only the equations (|60|l and (|62(l . Adding both 
equations in H(j2|) to the second equation in IjtiOfl gives (|65|) with p_i := 0. Hence, we can replace (|60ll - (|G2[1 with the 
linear system of equations consisting of the first equation in l|(j3D and the equations l|65l) , l(jtj|) . Now, by applying a 
similar reduction procedure as in the case k > 2, we obtain the polynomial 



P(l;X;(j,,u) =det I -1 A -1 I = A + 2(i = X + fi, 



1 
-1 A -1 
-2/i 2t- 2v A 

whose zero A = —fj, is the uniquely determined monodromy eigenvalue of A(^; /i, v) for each (fi, v) £ C 2 , ji 2 ^ v 2 . □ 
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Corollary 1 For a fixed half-integer k = k — i with a positive integer k, there exists a neighbourhood 11c C 2 of 
(0, 0) such that A(n; /z, v) has exactly 2k — I many monodromy eigenvalues Aq(k; [i,v), j = 1 — k, . . . , k — 1, for all 
(fx, v) £ it with /i 2 ^ v 2 . Moreover, Aq(k;/^, v) depends holomorphically on and lim^ lV )-t(o,a) Aq(k; Mi ^ = 3- 

In particular, monodromy eigenvalues and classical eigenvalues are distinct near (/i, v ) = (0,0). 

Remark 3 Monodromy eigenvalues also appear in the context of spheroidal wave equations. In EU Sec. 3.534] they 
are characterised by a similar property as given in Lemma\^ but they are not specified in detail. 

In view of Theorem [21 and Corollary ^ we could alternatively define the monodromy eigenvalues of A(k; /i, v) 
to be the zeros of the polynomial P(k; -\h,v) for each (li, v) £ C 2 (without the restriction /i 2 ^ v 2 ). Then the 
monodromy eigenvalues Aq(ac; 0, 0) — j, j = 1 — k, . . . , k — 1, fill in the gap of integers appearing in the spectrum 
of A(k; 0, 0). Moreover, P{\] A; [i, v) = A + \i is just the polynomial given by Lemma^Jin the case N = 1, and its 
zero Aq(j; /J,, v) = — A* satisfies the partial differential equation l|29|) for k = \. In the next section we prove that the 
monodromy eigenvalues of A{k; fi, v) are solutions of the PDE l|29|l for each half-integer n £ {i, |, |, . ..}. 



VI. MONODROMY PRESERVING DEFORMATIONS 

In 0, |10| and |ll| . Jimbo, Miwa & Ueno developed a general theory for monodromy preserving deformations of 
linear ordinary differential equations with rational coefficients. As a main result, they proved that the monodromy 
data (Stokes multipliers, connection matrices and exponents of formal monodromy) do not depend on the deformation 
parameters if and only if certain non-linear differential equations, the so-called deformation equations, are satisfied. 
This result, however, was proved under the restriction that the characteristic values at regular singular points do 
not differ by an integer. On the other hand, in the theory of special functions and in many physical applications 
the case where the characteristic values differ by an integer is of great significance. In this section we consider the 
isomonodromy problem for linear systems with two fixed regular singular points and coefficients which depend on one 
parameter t. Assuming that the characteristic values at the singular points are distinct and independent of t, we will 
show that certain components of the monodromy data are constant with respect to t if a deformation equation of the 
type |3, (1.18)] is satisfied. Since the monodromy components in question determine the existence of solutions of the 
form l|54(l and l|55|l . they are relevant to monodromy and classical eigenvalue problems. Applying the results to the 
system (II II) with an eigenvalue t of C as deformation parameter, it finally turns out that the deformation equation is 
in principle the characteristic equation of the partial differential equation 129|) . 

We start with some basic facts about parameter-dependent regular singular systems. At first, let us consider a 
family of (2 x 2) systems of differential equations 

^(x,t) = $(x,t)y(x,t), ( I ,f)e('B\{0})xS), (70) 

in an open disc 58 C C with centre that depends on a parameter t varying in some real or complex domain D. It 
is assumed that l|70|) has a regular singular point at for all t £ J). More precisely, we suppose that the coefficient 
matrix $ of i|70|) has the following properties: 

(a) &(x,t) = -^f(x,t), where : 55 x £> — > M2(C) is an analytical matrix function. 

(b) The eigenvalues a and j3 of ^(Q,t) are distinct and independent of t £ 53; moreover, Re a < Re/3. 

(c) There is an analytical function G : D — ► M2(C) such that G(t) is invertible and 

G(t) -1 *(0, t)G{t) = diag (a, (3) =: D, t£ D, 

Note that such a matrix function G always exists since the eigenvalues of ^(0, t) are distinct (see [3 Chap. VII, Sec. 
25, Theorem 25.1]). 

Lemma 6 If the conditions (a) - (c) are satisfied, then the system 11701 has a fundamental matrix of the form 

Y(x,t) = G{t)H{x,t)x D x J{t) 
where H : *B x D — ► M 2 (C) is analytic, H(0, t) = I for all t £®, and 

JW = U)o) « 71 » 

with some analytical function p : 5) — > C. Moreover, if (3 — a is not an integer, then p = 0. 
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Proof. If j3 — a is not an integer, the existence of such a fundamental matrix with J = is well known cf. 20]. 
Hence, we have to consider only the case that k := (3 — a is a positive integer. By the transformation 

y(x,t)=x a G(t)y (x,t), (72) 

the system (|7Uf) is equivalent to the differential equation 

x^(x,t) = * (x,t)y (x,t), (a; ) i)e(«8\{0})xD ) (73) 

where ^0(^7 1) '■— G(t)~ 1 ^(x 7 t)G(t) — a is an analytical matrix function, 

00 

$ {x,t) = J2 xn yo,n{ t )> (<M)g 05x25, 

?!=0 

with ^ofi(t) = diag (0, k) for all t G 25. Now, for j = 1, . . . , k — 1 we recursively apply the transformations 

= (^_ l(t) I) ^ 

where V'j-i denotes the (2, l)-coefHcient of the matrix ^j^i.i. At each step, yj(x, t) is a solution of a system 

x^(x,t) = 9 s (x,t)y 3 ix,t), (i,t)e(B\{0})xD, (75) 
where the coefficient matrix is analytic in 05 x 2), 



¥ i (a: > t) = £V*3.n(<). (M) 6 05 x 25, 

n=0 

with ^E , j,o(i) = diag (0, k — j) for all t € 2), and ^^(i), n = 1, ... , j — 1, are lower triangular matrix functions (that 
means, the (1, 2) component is identically zero). Finally, by the shearing transformation 

(l o\ 

Vk-iix-t) = Vk{x,t), (76) 

v° 

we obtain a differential equation 

=^k(x,t)y k (x,t), (x,t)€ (05\{0})x25, (77) 
where ^ : 05 x 2) — ► Ma(C) is an analytical matrix function, 

oo 

* fc (a?,t) = ^a; n * fc , n (t), (s.t) G 05 x 2), 

n=0 

satisfying 

**.o(*)=(^ t) °)=^), tes,, 

with some analytical function p : 2) — ► C. Note that p is just the (2, l)-component of ^k-i,i- Moreover, ^k,n(t), 
n = 0, . . . , k, are lower triangular matrices for alH G 2). Now, the system (|77|) has a fundamental matrix of the form 

Y(x,t) = H{x,t)x J{t) 

provided that H is a solution of the matrix differential equation 

OH 

x—(x,t)=y k (x,t)H(x,t)-H(x,t)J(t), (x,t) G 05 x 2), (78) 
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such that for each t G 35 the matrix H(x, t) is invertible for some, and hence all, x £ <8. Obviously, (|78() has a formal 
solution 

oo 

H(x, t) = ^T x n H n (t), (x, t) e <B x £, (79) 

n=0 

where H (t) = I and the coefficients H n , n > 0, are uniquely determined by the recurrence relation 

n-l 

(J{t) - n) H n (t) - H n {t)J{t) = -J2 *k,n-3(t)Hj(t) (80) 

3=0 

Following the proof of Theorem 5.3 in the book of Wasow |24|, it can be shown that the series (|79|l converges uniformly 
in every compact subset of 2$ x 2). Thus, a Weierstrass theorem implies that H is analytic in <B x 2), and therefore 
H is an actual solution of l|78|) . Further, since J(t) has the special form (|71|l and ^k,j(t), j — 0, . . . , fc, are lower 
triangular matrices, it is easy to verify that Hj(t) are lower triangular matrices for j = 0, . . . , k. Now, by combining 
the transformations 1(72(1 , (|74|l and ((76(1 , it follows that the differential equation 170(1 has a fundamental matrix of the 
form 

Y(x,t) =x a G(t) ( 1 ° ! ] H(x,t)x J W, (81) 



xq(x,t) x 

where q(x,t) is a polynomial in x of degree n — 1 with coefficients depending analytically on i, and H(x,t) is an 
analytical matrix function of the type 



H{x,t) 

satisfying hi%(0,t) — /i22(0, t) = 1. Now, if we define 
H(x,t) :-- 



hn(x,t) x k+1 h 12 (x,t) 
h 2 i(x 1 t) h 2 2{x,t) 



hn(x,t) xhi 2 (x,t) 
xq(x, t)hix(x, t) + x k+1 h 2 i{t) x 2 q{x, t)hi 2 (x, t) + h 22 (x, t) 



then H : <8 x D — > M 2 (C) is analytic, H (0, t) = I for all t <E D, and 

( / ° fc W(x ) t)=/f(a;,i)( J ° j. 
y xq(x, t) x k J \ x k J 

Hence, we can write the fundamental matrix 1(81(1 in the form Y(x, t) = G(t)H(x, t)x D x J ^ , where H has the properties 
stated in the Lemma. □ 

Now, we consider a family of (2 x 2) differential systems 

^(x,t)=${x,t)y(x,t), (x,t)e (<5\{0,l})xD, (82) 

in a domain (S, 2Jo U *Bi C © C C, with regular singular points at x = and x = 1 and a parameter t varying in some 
domain 5) C R or D C C. Further, we assume that the coefficient matrix <f> in ((82(1 has the form 

$(x, t) = , 1 ttfo i), (x, i) e (0 \ {0, 1}) x D, 
x{x — 1) 

where ^ : 25 x D — ► M2(C) is an analytical matrix function with the following properties: 
(I) The eigenvalues a, /? of ^(a,i) are distinct and independent of t € 5) and a G {0, 1}; in addition, Re a < Re/3. 
(II) There are analytical functions G a '■ 2) — ► M2(C), a € {0, 1}, such that G a (£) is invertible for alH 6 2) and 

G a (t)- 1 *(a,i)G a (<) = (-l) a diag(a,/?). 
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From Lemma it follows that the system (I82|l possesses a fundamental matrix of the form 

Y a (x,t) = G a (t)H a (x,t){x - a) D (x - a) J ^ (83) 

in the unit disc Q5 a C © with centre a £ {0,1}, where H a : <8 Q x 2) — ► M2(C) is an analytical matrix function 
satisfying H a (0, t) = I for all f £ 23, D = diag (a, /?), and 

J " <t)= (j*)o) (84 > 

with some analytical function p a : 2) — ► C. By analytic continuation along curves, we can assume that Y a is defined 
on the universal covering <R of the set © \ {0, 1}. Since Y a (x e 2i7r + a,t) = Y a (x + a, t)e 2i7rD [I + 2mJ a {t)\ for all 
(x, t) £ (Q3o \ {0}) x 2), the diagonal matrix D and the Jordan type matrix J a (t) represent the monodromy behaviour 
of Y a corresponding to a circuit around a £ {0, 1}. Moreover, as Yq and Y\ are both fundamental matrices of the same 
differential equation l|82(l . there exists an analytical matrix function Q : D — ► M2(C) such that Yq(x, t) = Y\(x, t)Q(t) 
for all (x, t) £ (<8\ {0, 1}) x 23, which is called the connection matrix for Yq and Y\. The next result gives a sufficient 
condition that certain components of the monodromy data J a and Q are constant in 2). For this reason, we establish 
in addition to (I) - (II) the following assumptions on the coefficient matrix $: 

(III) There exists an analytical function SI : © x 23 — ► Ma(C) such that 

^(x, t) + $(z, t)Q(x, t) = Q(x, t)$(x, t) + ^(x, t), (x, t) e (0 \ {0, 1}) x 2), (85) 

(IV) The matrix functions G a , a £ {0, 1}, satisfy the linear differential equations 

^(t) = il(a,t)G a (t), ten. (86) 

Theorem 4 If the conditions (I) - (IV) are satisfied, then 

dJ a _ dQ 2 i 



dt dt 



in 2), (87) 



where the Jordan matrices J a , a £ {0,1}, are given by (|84|l and Qyi : 23 — ► C denotes the (1,2) -component of the 
connection matrix Q for Yq and Y\ . 

Proof. Let 7 := /3 — a, and for fixed a £ {0, 1} we define 

QY 

Z a (x,t) := -gjr(x,t) ~ n(x,t)Y a (x,t), (x,t) £ m x 2D. 
From Q82fl and the deformation equation l|85|l it follows that 

dz a _ &Y a _ d n Ya _ ^ = 0(*n) _ *> _ 



dx dx dt dx dx dt dx 

and therefore Z a is a matrix solution of the differential equation (|82(l in 9L Hence, there exists an analytical function 
C a : 23 — ► M 2 (C) such that 

Z a (x,t) =Y a (x,t)C a {t), (x,t) £<Rx®. 
Now, by means of the differential equation l)86[l. we get 

^ = ^H a (x - a) D (x - a) J ° + G a ^(x - a) D (x - a) J " + log(x - a) G a H a (x - a) D ^{x - a) J « 
= (fi(a, • ) G a H a + Ga—Q^- + {x- a)~<log(x - a) G a H a -^j (x - a) D (x - a) J », 
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and since C a — Y a 1 Z a , it results that 

(x - a) D {x - a) J "C a (x - a)- ,a (x - a)-° 

= (x-a) D {x-a) J °Y- 1 (^-QY a ^ {x - a)~^(x - a)~ D 

f)M 3 7 

= H- l G-\n{a, • ) - fi) G a H a + H~ 1 -^- + (x - a) 7 log(x - a) 

= (x- a)F a + (x- ay \og(x -a)L J (88) 

with some analytical function F a : & x 2) — > M2(C). Further, by setting 

/ Cu(t) Ci 2 (i)\ 

aw -ia.( ( )a,(.)j- t€S ' 

(for clarity, we omit the index a in the entries of C a ), we have 
(x - a) D (x - a) Ja C a {x - a)- Ja {x - a y D 

/ 1 \ / C n C 12 \ / 1 \ 

y (x - a) 7 log(a; - a) p a (x-a) 1 J y C 2 i C22 J \-log(x-a)p a (x - a) -7 J 
_( Cn - \og(x - a)p a C 12 (x-a)^C 12 \ 
\ * C 22 + log(x - a) p a Ci2 



(89) 



Comparing (|88() to l|89|) . it follows that C12 = in £) since the function in l|88|l is bounded at x = a. This in turn 
implies Cn = C22 = as the diagonal entries in (|88|1 have a zero at x — a for all t € 35. Finally, we obtain that 

{x-a) D {x-a) J «C a {x-a)- J *{x-a)- D = I ° °] (90) 

V ; V ; y ' V ' \ [x-a)~<C 21 j y ' 

has no logarithmic singularity at x = a and therefore the last term l|88(l vanishes identically. Hence, ^f- = in D. 
Next, we prove that = 0. Since Y {x,t) = Yi(x,t)Q(t), it follows that 

dY °- dY i G + Y dQ foil 

Further, from Z a (x 1 t) = y a (x, t)C a (t) we get 

ay ay 

^^-^ri-yoCoyo- 1 , ^yf^-n^dYf 1 . (92) 



By means of <|91[1 and Y Q = Q 1 Y 1 , the first equation in (|92|l becomes 



^ y f 1 - n = y a q c q- 1 Yf 1 - y x ^ c^Yf 1 . (93) 



Now, (|9^|l and the second equation in (|9~2"|l imply 



Fi d yf 1 - Yi Q Co Q-%- 1 ~Yi^ Q- x Yy x 



and therefore 

dQ 



;)t =QC -C 1 Q. (94) 



Note that the matrix function C a has the form 



/ 0\ 



23 



Hence, if we set 



then l|94(l is equivalent to the system 



W)=: U(<>fe«> ' ,6S ' 



9 / On Q12 \ _ / c Qi2 \ 

dt \ Q21 Q22 J \ C0Q22 - ciQu C1Q12 J ' 

and we immediately obtain that d ^ 2 = in 55. □ 
In the following we apply the results of Lemma Eland Theorem 0] to a family of (2 x 2) differential systems 



dy t A 
d~x M = 



-B (t) + -1—B 1 (t) + C(t) 

X X — 1 



y(x,t), ( 3 ;,f)e(C\{0 ) l})x!D, (95) 



where f € J) with some domain 35 C K or 35 C C, and the coefficients Bo, B\, C : 35 — > M2(C) are supposed to be 
analytical functions. Further, we assume that the following conditions hold: 

(i) The eigenvalues a, (3 of B (t) are distinct and independent of t € 35. Moreover, they coincide with the eigenvalues 
of Bi(t), and Re a < Re/3. 

(ii) There are analytical functions G a : D — ► M2(C), a S {0, 1}, such that G a (t) is invertible and 

G (t)- l Bo{t)G (t) = -dtr^iWGift) = diag(a,/3) =: D, t€ 3). 

Let Y a , a £ {0, 1}, be fundamental matrices of (|95|l in the open disc Q5 a C C with centre a and radius 1 having the 
form iJSSJl, where H a : *B x D — ► M 2 (C) is analytical, H a (0,t) = I for all t G 35, and J a {t) is given by (JE1J with 
some analytical function p a : 3D — > C. Again, by analytic continuation, we assume that y a is defined on the universal 
covering $H of C \ {0, 1}, and we denote by Q : 3) — ► Ma(C) the connection matrix for Yq and Y\. In the sequel we 
are looking for conditions such that for fixed t € 35 the system l|95f) has one of the following properties: 

(P) There exists a fundamental matrix Y of the form 

Y(x) = [x(l - x)] a P(x)e c{t)x , xGC\{0,l}, (96) 

where P : C — > M 2 (C) is a polynomial matrix function. 
(H) There exists a nontrivial solution y of the form 

y(x) = [x(l-x)fh(x), xeC\{0,l}, (97) 

where h : C — ► C 2 is an entire vector function. 

Lemma 7 Suppose that the conditions (i) and (ii) are satisfied, and let t € 35 . TTien i/ie system l|95|) /ias i/ie property 
(P) z/ anrf onZy if (3 — a is an integer and pa(t) = pi(t) — 0, and it has the property (H) if and only if Q2i{t) = 0. 

Proof. By means of the transformation y(x) — x a (x — l) a y(x), the differential equation 195(1 is equivalent to the 
system 



d~x M = 



-Bo{t) + -^rBx{t)+C{t) 
X x — 1 



y(x,t), (x,t)€(C\{0,l})x 35, (98) 



where B>o(t) := B (t) — a and B>i(t) :— Bi(t) — a. Moreover, 



Y a (x,t)=G a (t)H a (x,t)(] ix _ a) p-a)( x - a y a(t) (") 

are fundamental matrices of l|98|) in a neighbourhood of a G {0, 1}. First, suppose that j3 — a is an integer and that 
Po{t) = Pi(t) — holds. In this case Jo(t) — Ji(t) — 0, and the system l|98|l has a holomorphic fundamental matrix 
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since (x — a)@ a is holomorphic and Y~o( ■ , i), Yx( ■ , t) contain no logarithmic terms. Moreover, as the coefficient matrix 
<!>( • , t) of l|95l) is a rational function which is bounded at infinity, the extension of Halphen's Theorem 0, Theorem 
2.4] implies that the system ()98(l has a fundamental matrix of the form Y(x) — R(x)e c ^ x with some rational (and 
hence polynomial) matrix function R : C — > M2 (C) . Conversely, if (|95|l has a fundamental matrix of the form l|96ll , 
then Yq( ■ , t) and Yi( • , t) are holomorphic matrix functions, which gives (3 — a £ Z and Jo(t) — J\{t) = 0. 
Next, let us assume that Q\2{t) = 0. If we define 

y{x) := Y (x, t) ( ° ^ = x^G Q (t)H Q (x, t) 

then y is a nontrivial solution of Ij95(l . and x~Py(x) is analytic at cc = 0. Since Fo(^,i) = Y\{x : t)Q{t) and Qi2(i) = 0, 
we obtain 

y{x) = Y 1 {x,t)Q{t) f ° J =(x-l) /3 G 1 (t)i/ 1 ( 2 ;,i) 

and therefore (1 — x)~^y(x) is analytic in a neighbourhood of ir = 1. Now, by the existence and uniqueness theorem, 
h(x) := — x)]~' 3 y(a;) can be extended to an entire vector function. Conversely, suppose that l|95l) has a nontrivial 
solution of the form (|97|l . Then 

y(a;) = F (z,t) ( ° J = Y x {x,t) ^ ° ^ 
with some constants Co, ci € C \ {0}. Since Yo(x,t) — Y\(x,t)Q(t), it follows that 

Q(t) 






which gives Q2\{t) =0. □ 
Now, in addition to (i) - (ii), we assume that the coefficients of (|95|l satisfy the following conditions: 

(iii) There exists an analytical function f2 : C X 3) — ► Ma(C) such that the deformation equation (|85|) holds in 
(C \ {0, 1}) x 2), where $ is given by 

$(z,i) := - B (t) + Bt(t) + C(t), (i,i)£(C\{0,l})x2. 
x x — 1 

(iv) The matrix functions G a , a £ {0, 1}, satisfy the differential equations 

^(t) = n(a,t)G a (t), tea). 

The next result is an immediate consequence of Theorem 0] and Lemma 

Corollary 2 Suppose that the conditions (i) - (iv) are satisfied. If (P) holds for one to e 55, t/ien <|95|) /ias t/ie 
property (P) /or all t £ 2), and z/ (H) Zio/ds /or one to G 2), i/ien (|95|) /ias the property (H) /or all t £ D. 

Finally, we apply the results of this section to prove that the classical as well as the monodromy eigenvalues of the 
Chandrasekhar-Page angular equation in dependence of (/z, v) € R 2 are (locally) solutions of the partial differential 
equation jSSjl. 

Lemma 8 Let k £ [5,00) and a e {— 1, +1} be fixed. Moreover, suppose that the functions v : 2) — ► K \ {0} and 
w : D — > R are solutions of the system (|50|l - <|51[1 on some interval 2) C (0, 00). Finally, let 

4) > (ioo) 

and to £ 2). If w(to) is an eigenvalue of A(k; /j,(to),v(t )), then w(t) is an eigenvalue of A(k; fj,(t) , u(t)) for each 
t £'D. Furthermore, if k is a half-integer and to (to) is a monodromy eigenvalue of A («;; /i(to), v(to)), then w(t) is a 
monodromy eigenvalue of A(n; jj,(t),v(t)) for each t £ 2). 
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Proof. In terms of (|100fl and A(i) :— w(t), the coefficient matrices (|12|) of the system take the form 



- + - 




1 o t+i v H w(t) + =w) - -s - * y 

-(«(<) 2 +i) -(^) 2 -i)\ 

( W (t) 2 -i) («(i) 2 +i)/' 

where the condition (i) is satisfied with a = — 5 — | and /? = § + z- Now, if we define the matrix function 

l f( v (tf-a)(i-x) (v(tf + a){l-x)\ 
n(x, t) := — — — I " , (i,i)eCxD, 

"(*) V {v(tf + a)x (v(tf - a) (x - ±) ) 

then, by a straightforward calculation using the characteristic equations (|50|l - J5TJ, it follows that the deformation 
equation in (iii) holds. Finally, by setting 




(«(*) + w{t) 
( K+ i) e -*(t) 



-0(i) 



, Gi(i) := KG (t)K, 



where if is the matrix l|20(l and 



J to 



2«(r) 



the conditions (ii) and (iv) are satisfied. Since a point w(t) is a monodromy eigenvalue of A (k; ^(i)) if and only 
if has the property (P), and w(t) is a classical eigenvalue of A (k; n(t), v(t)) if and only if has the property 
(H), the assertion follows from Corollary [21 □ 

Theorem 5 For a fixed K = k — \ with a positive integer k, let (0, 0) £ 6 C I 2 be a simply connected domain such 
that for each (/i, u) € 6 all monodromy eigenvalues Aq(k; fi, v), j = 1 — k, . . . , k — 1, of A(k; fi, v) are simple zeros of 
the polynomial P{k; ■ , [i, v) given by Theorem^ Then each function A = Aj, j = 1 — k, . . . , k — 1, satisfies the partial 
differential equation (|29|l in 6. 

Proof. Let j € {1 — k, . . . , k — 1} be fixed. The monodromy eigenvalues of A(k; /j, v) are exactly the zeros of the 
polynomial P(k; ■ ; [i, v), and since all zeros of P{k; •;h,v) are simple, the implicit function theorem implies that 
Aq(k; /x, v) depends analytically on (fi, v) in &. In order to show that the function A = Aq satisfies the PDE l|29|) . we 
make use of the unique continuation property of analytical functions. That means, it suffices to prove that (|29|l holds 
for A = Aq in a neighbourhood of some point (//, v) = (r, 0) € 6, r > 0. Now, in view of the coordinate transformation 
(|100[1 . we have to verify that the function Aq (k; fj,(t, v), v(t 1 v)) is a solution of the partial differential equation 14911 in 
a neighbourhood of the point (t, v) = (r, 1). To this end, let us consider the characteristic equations of l|49|l 

dv 2v(t,u)w(t,u) dw ( 1 \ t ( 2 1 



-(t,u) = , -5r(*,u) = -K V(t,u) + — r - - v(t,U 



dt y ' ' t ' dt y ' ; V ' u 0> u ) / 2 V u ) 2 

together with the initial values 

v(t, u) = u, w(t, u) = Aq (k; jLt(r, u), i/(t, it)) , 

which depend analytically on the parameter u G (0, oo). The solutions v(t, u) and u) of this initial value problem 
are analytical functions in a neighbourhood of (t, 1), and since ^(r, u) — 1, they form locally an integral surface 
for the PDE (|49J) (compare [l2l Chap. 1, Sec. 5]). More precisely, there exists an analytical function U defined on 
a neighbourhood 03 of (t, v) = (t, 1) such that U(t,v) = v, and W(t,v) :— w (t,U(t,v)) is a solution of Ij49(l in QJ. 
Now, Lemma [S] implies that W(t,v) is a monodromy eigenvalue of A (k; /i(t, v), v(t, v)) for all (t, v) G 2J, and since 
W^r, w) = Aq (k; /x(r, u), ^(r, u)), it follows that W(t, v) = Aq (k; /i(t, v), v(t, v)) holds identically on 33. This completes 
the proof of the Theorem. □ 
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In a similar way we can apply Lemma[S]to prove that for fixed k S (0, oo) the zeros of the function A i — > A(k; A, /j,, v) 
defined in Section II and therefore the eigenvalues of A(k; /z, v) satisfy the partial differential equation 129( 1 . This 
alternative proof of Theorem ^ is based on monodromy preserving deformation - a general technique, which should 
be applicable to other eigenvalue problems as well. Potential candidates and associated fi-matrices for solving the 
deformation equations can be found in H3, Appendix C]. 

Finally, as a consequence of Theorem the zeros of the polynomial P(k; ■ ; fi, v) given by Theorem |21 satisfy the 
PDE l(29|) and do not coincide with any eigenvalue of A(k; fi, v) in a neighbourhood of (fx, v) = (0,0). Moreover (see 
the proof of Lemma^J), P(k; ■ ; fi,v) gives rise to a special integral of polynomial type for the Painleve III <(52|) . Now, 
the results of Mansfield & Webster in 0, Section 2] suggest that these special integrals are unique in some sense, 
which in turn implies that classical eigenvalues of the Chandrasekhar-Page angular equation are not algebraic. 

Appendix 

A. Eigenvalues and eigenfunctions in the case fi = v = 

For fixed n £ [i, oo), a point A £ C is an eigenvalue of A(k; 0, 0) if and only if the system 10 with (p, v) = (0, 0) 
has a nontrivial solution S(9) satisfying 



/•7T 

/ \S(9)\ 2 d6 < oo. 
Jo 



(101) 



Introducing the functions u, v : (—1,1) > C by 



'tan | w(cos6 

S{6) =: sin K +t I \ ' I , (102) 



' cot | w(cos ( 



then JSJ) with (/i, v) — (0, 0) is transformed into 



(1 - x)u'(x) = + u(x) + Xv(x), (1 + x)v'(x) = -Xu(x) ~ {^ K +\) riJ ' 1 " |: " !) 

and the normalisation condition i|101fl is equivalent to 

l pi 



u{x) z {l-x) K+ i{l + x) K -i Ax < oo, / v(xy(l-x) K -i(l + x) K+ 2 dx < oo. (104) 
l J-i 

If A = 0, then the differential equations (|103fl imply u(x) — Ci(l — x)~ K ~^ and v(x) = C2(l + x)~ K ~^ with some 
constants C\, c% € C, and from the condition (|104fl it follows that ci = c 2 = 0. Hence, A = is not an eigenvalue of 
A(k; 0, 0), and we assume in what follows that A ^ 0. In this case, the second equation in (|103|) gives 

u{x) = -^-^v'(x)-^^-v(x), (105) 
A A 



and for v we obtain the second order differential equation 



(1 - x 2 )v"(x) + [1 - 2(k + l)x] v'(x) + 



V 2 



v(x) = 0. 



If we set ^,/9:=K+i and A := A — k — i, this differential equation becomes 

(1 - x 2 )v"{x) + [f3-a-{a + f3 + 2)x) v'[x) + A(A + a + (3 + l)v(x) = 0, (106) 
and the second condition in (1 104(1 takes the form 

v(x) 2 (l-x) a (l + x) p dx < oo. (107) 
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Note that (|106fl and i|107|l is the eigenvalue problem associated to the Jacobi polynomials. More precisely, the solutions 
of the differential equation IjlObj) which are square integrable with respect to the weight function (1 — x) a (l + x) 13 

are constant multiples of the Jacobi Polynomials Pn with some non-negative integer n, and the corresponding 
eigenvalues are determined by the equation A 2 — (k + i) = n(n + a + (3 + 1), i.e., A^ = ± (k + \ + n). Now, if 
we define v(x) := -Pn a,l3 \x), x G (—1, 1), then l|105[l yields 

A± u(x) = (1 + X)±PM +fU**f» = " + + + x)PtX U3+1) +PPt' l3) 

+ f3Pt?' 0) 

= (a + l3 + n + l)p( a+1 <V ~(a + n+ l)P 7 [ aJi) = (/? + n)P 7 [ a+1 ' f3 -^ = |A± \p(«+W-i) 

where we applied the differentiation formulas and contiguous relations for Jacobi polynomials (see [Tol Section 5.2]). 
Hence, u{x) = ±P I ( " +1 ' /3 ~ 1) (x), x £ (-1, 1), and since u satisfies the first condition in (|f 04fl . the numbers A„ are in 
fact eigenvalues of A(k; 0, 0). Moreover, the corresponding eigenfunctions are constant multiples of 

sin«0 v 2 , ee(0,7r), 

\ -yW|P r ( r 5A+5) (cos^) J 

which form a complete orthogonal set in Jzf 2 ((0,7r),C 2 ). In particular, the spectrum of A(k;0, 0) is given by {A^ : 
? i = 0,l,2,...}. 



a + p + n + 1 
ni.+ R + In + 1 



{P + n)Pt\^ ] +nP^ 



B. A Numerical Example 

As a numerical example, we have computed the coefficients c m , n of the power series expansion 1)38(1 up to and 
including m + n — 8 for n — | and j = 1 using the recurrence relation given in Section III. The coefficients have 
been rounded to six significant figures and listed in the table below. It should be noted that they are to some extent 
different from the coefficients displayed in [2^, Table I]. Evaluating the power series expansion l|38l) at a = 0.01 and 
(3 = 0.02, i.e., (/!, v) = (0.005,0.015), yields Aj_ = 1.01167 as a numerical approximation for the eigenvalue Ai, and 
this result coincides with the value given in |23j, Table II]. For a second pair of parameters (a, (3) = (0.5,1.0), i.e., 
= (0.25,0.75), we obtain A\ = 1.59745, which differs slightly from the value Ai = 1.59764 listed in JH, Table 
II]. In order to test the reliability of our numerical result, we can use the statement of Lemma |3| That means, we 
approximate 8(A) defined in H25|) by the second component 0„(A) of d n (X) for n = 8, and we compare 0g(Ai) and 
6 8 (Ai) with the theoretical result 9(Ai) = 0. As 6 8 (Ai) = 3.60882e - 05 and 6 8 (Ai) = -2.51164e - 04, our result 
seems to be more trustworthy. Finally, let (/j,, v) — (0.02,0.1). The coefficients of the polynomial 8 are given in 
Table ITT1 For these parameters, our power scries approximation gives X\ = 1.07379 which differs significantly from 
the value Ai = 1.06104 given by Chakrabarti (see 0, Table 1]). Despite his claiming of an accuracy of six decimals, 
the evaluation of 8 at the eigenvalues in question gives 9 8 (Ai) = 5.68899e — 12 and 6 8 (Ai) = 1.52770e — 02 in 
favour of our result. Thus, Chakrabarti's calculations should be taken with some caution. 







1 



6 



7 



8 



n = 
1 
2 
3 
4 
5 
6 
7 



1.00000e+00 
1.66667e-01 
7.40741e-02 
-8.23045e-03 
-9 . 14495e-04 
.08053e-04 
.38702e-05 
.63435e-05 
,10856e-06 



5.00000e-01 
0.00000 

-1.48148e-02 
3.29218e-03 
5.48697e-04 

-4.06442e-04 
3.38702e-05 
3.16122e-05 



0.00000 
0.00000 
0.00000 
-4.70312e-04 
-1.22281e-04 
1.41790e-04 
-1.63351e-05 



0.00000 
0.00000 
0.00000 
0.00000 
35868e-05 
67091e-05 



0.00000 
0.00000 
0.00000 
0.00000 
0.00000 



0.00000 
0.00000 
0.00000 
0.00000 



0.00000 
0.00000 
0.00000 



0.00000 
0.00000 



0.00000 



TABLE I: The coefficients c m , n ,0<m + n<8, of the power series expansion 1381 in the case k = i and j = 1. 
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n = 


1 


22151e+00 


n = 9 


4 


91151e- 


-06 


1 


1 


44347e-02 


10 


-4 


22048e 


-04 


2 


-1 


70525e+00 


11 


-9 


46610e 


-08 


3 


-7 


92297e-03 


12 


1 


02643e- 


-05 


4 


6 


70 AAA— r\ H 

72114e-01 


1 9 

lo 


6 


88933e- 


-10 


5 


1 


46003e-03 


14 


-1 


26470e 


-07 


6 


-1 


12351e-01 


15 


-1 


OOOOOe 


-26 


7 


-1 


21028e-04 


16 


6 


15119e- 


-10 


8 


9 


39664e-03 











TABLE II: The coefficients S n of the polynomial 6 8 (A) = J2n=o a « A " for K = §> A* = 0.02, i/ = 0.1. 



C. Eigenfunctions in the case 7^ |f| 

Eliminating the second component of y in the system (lllf) . we get a linear second-order differential equation for 
the first component y\ given by 

^ ( " )+ U-^J^ (s)+ ( v TO + - + ^ + ^ + (^IF + ^J yi(a;) = (1 ° 8) 



with 



and 



b:=?-z , r := 4 (/x 2 - z/ 2 ) , n := A 2 - 2 a 2 + 2 v + a - /i 2 - 4 a v + 



a := - + - 



2/i ^ / ' - f ix— \ 2 4 



2 4 a// 2 2(^ 2 + 2a/ f 2 - 1 /A 2 ) 

r 2 := -a , r 3 := — — +2i/-ti, r 4 := a(l - a), t 5 := -5 . 

/i z — \ z X z — \i A 

Now by means of the transformation 



yi (x):=x a (x-l) a ip(x)e 2tx , t = ±V^ 2 - A* 2 - 
we find that VK^) satisfies the generalised Heun equation 

d 2 ip(x) | / 1-Mn | 1 — A*i | 1-/^2 | #(>) | A) + As + P2X 2 ^ _ Q , 1Qg . 
dx 2 \ ^ x— 1 x — & J dx x(x — l)(x — b) 

where 



and 



^ = -2a, /Lt x =l-2a, ^ 2 = 2, /3 2 :=8ai, 



, o , ,2a tt(6 — 4) 2a ab 



/?! = /i 2 - A 2 - 2 i [b + 2 a(l + 2 6)] + 2 a(a> - 4) + 2 z/(2 a - b) - 
f3 a = b (A 2 - ii 2 ) + b[2{v + t) - Aa{v - t) - Aa 2 } + a - 



A + /! /! — A 

2 /i a 6 



A — /! 

We observe that 0, 4 and b are simple singularities with characteristic exponents (0, /io), (0, fii) and (0, ^,2) respectively, 
while 00 is (at most) an irregular singularity of rank 4. To stress the importance of equation H409(l . it is sufficient to 
remark that it contains the ellipsoidal wave equation as well as Heun's equation and thus the Mathieu, spheroidal, 
Lame, Whittaker-Hill and Ince equations as special cases. 
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